Universal probes of two-dimensional topological insulators: Dislocation and vr-flux 



Vladimir Juricic,^ Andrej Mesaros,^^ ^ Robert- Jan Slager,^ and Jan Zaanen^ 

^ Instituut-Lorentz for Theoretical Physics, Universiteit Leiden, 
P.O. Box 9506, 2300 RA Leiden, The Netherlands 
'^Department of Physics, Boston College, Chestnut Hill, MA 024 67, USA 

We show that the 7r-flux and the dislocation represent topological observables that probe two- 
dimensional topological order through binding of the zero-energy modes. We analytically demon- 
strate that TT-flux hosts a Kramers pair of zero modes in the topological F (Berry phase skyrmion 
at the zero momentum) and M (Berry phase skyrmion at a finite momentum) phases of the M — B 
model introduced for the HgTe quantum spin Hall insulator. Furthermore, we analytically show that 
the dislocation acts as a vr-flux, but only so in the M phase. Our numerical analysis confirms this 
through a Kramers pair of zero modes bound to a dislocation appearing in the M phase only, and 
further demonstrates the robustness of the modes to disorder and the Rashba coupling. Finally, we 
conjecture that by studying the zero modes bound to dislocations all translationally distinguishable 
two-dimensional topological band insulators can be classified. 



The topological band insulators (TBIs) in two and 
three dimensions have attracted a great interest in both 
theoretical and experimental condensed matter physics [1, 
[5]. An extensive classification of topological phases in 
free fermion systems with a bulk gap, based on time- 
reversal symmetry (TRS) and particle-hole symmetry 
(PHS), has been provided [3], and it defines the so-called 
tenfold way. This however relies on the spatial continuum 
limit while TBIs actually need a crystal lattice which, 
in turn, breaks the translational symmetry. Therefore, 
the question arises whether the crystal lattice may give 
rise to an additional subclass of topological phases. The 
simple M — B model introduced for the two-dimensional 
(2D) HgTe quantum spin Hall insulatorf^ gives away a 
generic wisdom in this regard. Depending on its parame- 
ters, this model describes topological phases which are in 
a thermodynamic sense distinguishable: their topological 
nature is characterized by a Berry phase skyrmion lattice 
(SL) in the extended Brillouin zone (BZ), where the sites 
of this lattice coincide with the reciprocal lattice vectors 
('T-phase") or with the TRS points (tt, tt) ("M-phase"). 

The question arises how to distiguish these phases by 
a topological observable. Early on it was observed in 
numerical calculations that a TRS localized magnetic tt 
flux binds zero- modes in the r-phase[5]-[7]- It was also dis- 
covered that screw dislocation lines in three-dimensional 
TBIs carry zero mode structure's]. Here we will demon- 
strate the special status of such topological defects as 
the universal bulk probes of the electronic topological 
order. We will present an analytical description of the 
zero modes bound to tt flux, demonstrating that these 
are closely related, but yet different from the well-known 
Jackiw- Rossi solutions [5]. These modes also generalize to 
two-dimensions the effect of one-dimensional spin-charge 
separation [in]. By performing numerical computations, 
we will show that the dislocation binds zero modes only in 
the topologically nontrivial M-phase. Subsequently, we 
will demonstrate that, modulo a change of basis, these 
dislocation zero modes are of the same kind as the tt flux 



ones, but appear only in the case of non-F-type TBIs. We 
conjecture that by studying the presence of zero modes 
associated with dislocations all possible "translationally 
distinguishable" TBIs can be classified at least in the 2D 
case. 

Consider a tight-binding model proposed to describe 
HgTe quantum wells [4] 
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where — {u-ff,v^,u^,v^) = (^'^,'1'^), while u and v 
represent two low-energy orbitals. El and HI in the case 
of the HgTe system. The upper and the lower blocks in 
the Hamiltonian are related by the time-reversal symme- 
try, and i?(k) acting in the orbital space has the form 

F(k) = T^d^(k), (2) 

where t^, /i = 1,2,3, are the Pauli matrices, di.2(k) — 
sin k^ y, and = A/— 2i?(2— cos fca;— cos ky), while lattice 
constant a — 1, h = c — e — 1, and summation over 
the repeated indices are assumed hereafter. The Dirac 
mass M and the Schrodinger mass B are related to the 
material parameters, while the form of the Hamiltonian 
([l]) is set by the symmetries of the system [TT]. Since 
the above Hamiltonian has the spectrum i?(k) = ^Jd^d^ 
doubly degenerate in the spin space, the band gap opens 
at the F-point for < M /B < 4 and the system is in 
a topologically nontrivial F-phase. For 4 < AI/B < 8, 
the system is still topologically nontrivial, but with the 
bandgap opening at k = (tt, 0) and (0, tt) (M-phase). On 
the other hand, for M/B < and M/B > 8 the system 
is topologically trivial with the bandgap located at the F 
and the M point, respectively. 

In the topologically nontrivial F-phase, the band- 
structure vector field d(k) = d(k)/|d(k)| [Eq. ([2|] forms 
a skyrmion centered at the F-point in the BZ, Fig. [l|^a) 
for AI/B ~ 1, with the corresponding skyrmion density 
s(k) = d(k) • (afc^d(k) X dk^d(k)) shown in Fig. Qb) 



[s(k) tracks the position of minimal bandgap in the BZ, 
coinciding with it where d(k)||(9j.^d(k) x i9j.^d(k)]. In 
the 2D extended BZ, this skyrmion structure forms a 
lattice which respects point group symmetry of the orig- 
inal square lattice. Furthermore, in the M-phase, the 
skyrmion is centered at the M point in the BZ, Fig. 
[ija) for M/B = 5. The position of the corresponding 
skyrmion lattice relative to the extended BZ is therefore 
different than in the F-phase, but the skyrmion lattice 
still respects the point group symmetry of the square 
lattice. On the other hand, in the topologically trivial 
phase, the vector field d forms no skyrmion in the BZ, 
as shown in Fig. [ija) for M/B = -1 and M/B = 9, 
consistent with the vanishing of topologically invariant 
spin Hall conductance crf^ = (47r)^^ /^^ d^k s(k). The 
position of the skyrmion lattice relative to the extended 
BZ thus encodes translationally active topological order 
which, as we show below, is probed by the lattice dislo- 
cations. 

Let us first numerically demonstrate that lattice dislo- 
cations bind zero modes as robust topological phenomena 
in the Af -phase of the M — B model. Analysis in Ref. [12] 
predicts that PHS (which our model respects), topologi- 
cally protects localized zero modes on dislocations when 
the system has K^dge = edge modes which our model 
has only in the M phase. We indeed find in the M-phase 
a Kramers pair of zero modes bound to a dislocation 
which results from the fact that, as we show below, the 
dislocation acts as a tt-Aux in this phase. However, we 
here show the robustness of these modes in the bulk gap 
when we introduce the random chemical potential, and 
thus break PHS. To this end we performed numerical 
analysis of the tight-binding M — B model in the real 
space 



Htb — 




T5 + z^(]H-T3)e,.(<Tx5) 



where — (s^(R),p^(R), s^(R),p|(R)) annihilates 
the \E1, ~ \s) type, and \H1, |) ~ \px + ipy) type or- 
bitals at site R and nearest neighbors 5 € {e^,, Ej,}; Pauli 

matrices cr mix spin. We set Ts.^-f = ^ , Tg^n — 

Tg^^^, with tx=t'^^ -i, ty = ~t'y = -1, A^/p = ±B+D 
and on-site energies e = [(C — 4£')to + {M — 4_B)r3] ® ctq. 
This reproduces Eq. ^ when C ^ D = {), implying 
vanishing chemical potential. The Rq term is the nearest 
neighbor Rashba spin-orbit coupling [13) which is induced 
by broken 2 — )• — z reflection symmetry of the quantum 
well, and should be relevant in case of tunneling mea- 
surements on thin wells. 

Our numerical analysis of Htb pertains to various sys- 
tem shapes and sizes, with varying disorder strengths 
given by multiplication of the parameters for each R, 5 



by Gaussian random variables of width w, while preserv- 
ing TRS. Fig.[IJc) demonstrates the spectrum and wave- 
function at w = 10% disorder. Localization (even by 
weak disorder) decouples dislocation states from edges 
and possible edge roughness effects. 

In the presence of the Rashba coupling (i?o ^ 0, but 
not large enough to close the topological bulk gap 14 ), 
the spins are mixed, but the Kramers pairs remain local- 
ized (Fig.[T]i). Fig. [2] demonstrates the robustness of dis- 
location modes within the topological bulkgap, through 
the Rashba coupling perturbed and disorder averaged 
density of states (DOS) of Htb in a periodic lattice, con- 
trasted between the F and the M phase. 

Let us now analytically show, using the elastic contin- 
uum theory, that the effect of a lattice dislocation in the 
M-phase is to effectively introduce a magnetic tt flux. We 
consider a dislocation with Burgers vector b, and expand 
the Hamiltonian ([2| around the M-point in the BZ. As 
a next step, a dislocation introduces an elastic deforma- 
tion of the medium described by the distortion {ei(r)} of 
the (global) Cartesian basis {e^}, i = x,y, in the tangent 
space at the point r|15|, the momentum in the vicinity 
of the M-point reads 

k, = Ei • (kM - q) = {ei + £i) • O^M - q), (4) 

where k^f = (7r,7r)/a, q is the momentum of the low- 
energy excitations, |q| ^ |k7\/|, and we have restored the 
lattice constant a. The corresponding continuum Hamil- 
tonian after this coarse graining|16j is 

H,s{k,A)^n{h+A,) + [M-B{k + A)^]T3, (5) 

with redefinition q k, M = M — 8B, B = —B, and 
Ai = —Si ■ kif. The form of the distortion £i is deter- 
mined from the dual basis in the tangent space at the 
point r which in the case of a dislocation with b = ae^ 
is = (1 - 2^)e- + ^ey,-Ey = e^^fTT]. Using that 
Ei • E^ — Sf, we obtain the distortion field to the leading 



(3) order in a/r, 



Finally, this 



form of the distortion yields the vector potential 



2r2 



(6) 



in Eq. ([5|, and therefore the dislocation in the M-phase 
acts as a magnetic 7r-flux. On the other hand, in the F- 
phase, the continuum Hamiltonian has the generic form 
|5|, with M = M and B = B, but the action of the dislo- 
cation is trivial, since the bandgap is at zero momentum, 
and thus A = 0. 

The last operation is to show the origin of the zero- 
modes bound to the 7r-flux/dislocation. This involves 
a generalization of the Jackiw- Rossi (JR) mechanism[9]. 
We flnd that instead of the vortex in the mass required 
for the JR solutions, quite similar zero modes are formed 
by fermions with a momentum-dependent mass in the 
background of a 7r-flux. In the F phase, we express the 
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FIG. 1. (color online) Model of Eq. ([l]) in the BZ: (a) Band-structure d(k). (b) Skyrmion density s(k). Mid-bulkgap localized 
dislocation states in 33x30 unit-cell M — B tight-binding lattice with disorder. The Kramers degenerate pair states are omitted, 
(c) Dislocation in the center. Offset disks represent the amplitude of s tiP t states, and the color their phase, (d) Total 
wavefunction amplitude in a periodic system (necessitating two dislocations), with Rashba coupling (-Ro) mixing spins. 




FIG. 2. Comparison of F (a-d) and M (e-h) phases. The density of states of 21x18 lattice (100 disorder realizations averages), 
with C = 0-2\t\, D = 0.3\t\ setting the chemical potential, and Rq the Rashba coupling, (a) and (e) In absence of dislocation. 
(b,f) Robust midgap dislocation modes are present only in M-phase; (c,g) Same is true upon spin mixing through Ro- (d,h) 
Strong Rashba coupling closes the topological bulkgap. 



Hamiltonian ^ in polar coordinates {r,ip), with M — 
M, B = B, the vector potential as in Eq. ([g]), and use 

the ansatz for the spin up zero-energy state 



equation: 



-h {2MB ~l)Oi_i 



B^Of_i 



w/_i(r) = 0, (8) 



ui-i{r) 



(7) 



with ^ G Z as the angular momentum quantum num- 
ber and the spin index for the spinor suppressed. The 
function u is then found to be solution of the following 



where the operator Oi = + r ^dr — r . The func- 
tion vi (r) obeys an equation of the same form as ^ with 
l->-l + l. 

This result also follows by noting that if the spinor in 
Eq. ([t]) is an eigenstate with zero eigenvalue of Hcs (k. A) , 
then it is also an eigenstate with the same eigenvalue of 



the square of i?cff(k, A). One then obtains 

iJeff(k, A)2 = B^{k^f + (1 - 2MB)k^ + M^, (9) 

with k = k + A, and the operator k^ after acting on the 
angular part of the upper component of the spinor Q 
yields Eq. ([S]). From Eq. ^ we conclude that the func- 
tion w/_i(r) is an eigenfunction of the operator Oi_i/2 
with a positive eigenvalue 



0/_iu;-i(r) = A^M/_i(r), 



(10) 



since the operator k^ when acting on a function with 
the angular momentum I is equal to — C/+1/2, and the 
eigenstates of the operator k^ with a negative eigenvalue 
are localized. Eqs. Q and ( 10 ) then imply 



Ah 



1 ± Vl - 4MB 
2B ' 



(11) 



and the function ui (r) is a linear combination of the mod- 
ified Bessel functions of the first and the second kind, 
/;_i/2(Ar) and Ki_i/2{^i'), respectively. Their asymp- 
totic behavior at the origin and at infinity then implies 
that the only square-integrable solutions are in the zero 
angular-momentum channel. Note that the localization 
lengths \± coincide with the penetration depths asso- 
ciated with the edge modes. We should now distin- 
guish two regimes of parameters, < MB < 1/4 and 
MB > 1/4, for which the argument of the square-root in 
Eq.( |ll[ ) is positive and negative, respectively. 

For < MB < 1/4, A± are purely real, and using Eq. 
([8|) we obtain zero-energy solution in the form 



vl/(r)-vl/^(r) + *_(r). 



where 
*±(r) 



C^e"^^-"-) 



±i 



(12) 



(13) 



with complex constants. For both M and B positive 
(negative), A± > (A± < 0), and therefore only the 
solution (5*-) is normalizable. For MB > 1/4, the 
existence of the zero-modes can be similarly shown. 



The zero-energy modes ( 13 1 form an overcomplete ba- 



sis in the zero angular-momentum channel as a conse- 
quence of the singularity of the vortex potential. A par- 
ticular regularization provided by considering the vortex 
with the flux concentrated in a thin annulus ensures Her- 
mitianity of the Hamiltonian [T51[TH]. This regularization 
yields, up to normalization, 



vl'(r) 



(14) 



Notice that this zero-energy state is regular at the origin 
which is a consequence of the regularity at the origin of 



the solutions in the absence of the vortex. Other regu- 
larizations, representing different microscopic conditions 
at the vortex, can lead to different behavior [20] . 

This result for the bound states to the vortex in the 
quantum spin Hall state should be compared with the 
one for a trivial insulator when MB < 0, and to be con- 
crete we consider B < 0. Then, A+ < and A_ > 0, 
and normalizable solutions are given by Eqs. ( 12 ) and 
( 13 1 with Ci ~ C2 = 0. However, since the upper and 
the lower component of the spinor diverge at the origin 
with opposite signs, the condition of regularity cannot be 
satisfied simultaneously for both, and therefore no zero- 
energy mode exists in the topologically trivial phase. No- 
tice that the condition of regularity at the origin is anal- 
ogous to the boundary condition for a TI interfaced with 
the trivial vacuum|llj , while the flux provides topological 
frustration to the electrons in the bulk. 

These results also pertain to the 7r-flux in the topo- 
logically nontrivial M-phase M/B > in Eq. ([s]), now 
involving the coarse grained states around the M point, 
where we also find zero-modes bound to 7r-flux. Since the 
dislocation in the Af-phase, as we have shown, effectively 
acts as a tt-Aux, it also binds a pair of zero modes in this 
phase, and in fact, this explains our numerical resuls. On 
the other hand, in the trivial phase M/B < 0-neither the 
TT-flux nor the dislocation binds zero modes. 

Depending on their occupation, the modes bound to 
the topological defects carry nontrivial charge or spin 
quantum number. The spin-charge separation, character- 
istic for one-dimensional systems [10], thus appears also in 
a two-dimensional system, and is tied to a topologically 
non-trivial nature of the quantum spin Hall state [S] [7]. 

These findings are obviously consequential for experi- 
ment. Dislocations are ubiquitous in any real crystal, for 
instance in the form of small angle grain boundaries. We 
predict that their cores should carry zero modes, which 
should be easy to detect with scanning tunneling spec- 
troscopy. The experimental challenge just lies in the re- 
alization of non-F TIs that are also easily accessible to 
spectroscopic measurements. 

Based on these specific results we conjecture the follow- 
ing general principle. Besides the "translationally triv- 
ial" F-type TBIs, completely classified by the tenfold way, 
there is a further subclassification in terms of "transla- 
tionally active" TBIs which are in 2D characterized by 
the locus of the Berry phase skyrmion lattice in the ex- 
tended BZ, involving high symmetry points respecting 
time-reversal invariance and point group symmetries (like 
the M phase). The 7r-fluxes and dislocations act as the 
topological bulk probes and the presence or absence of 
zero modes can be used as a classification tool for the 
translational side of topological order in TBIs. For in- 
stance, the Haldane phase on the honeycomb lattice [21] 
has skyrmions at the valleys and accordingly zero modes 
bound to dislocations. At present we are testing this 
hypothesis for all wallpaper groups in 2D. 
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